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Abstract 

Using Nambu^spin space Keldysh Green's function approach, we present a nonequihbrium 
charge and spin pumping theory of a quantum dot in the mico-cavity coupled to two superconduct- 
ing leads. It is found that the charge currents include two parts: The dissipationless supercurrent 
standing for the transfer of coherent Cooper pairs and the pumped quasi-particle current. The 
supercurrent exhibits a dynamic — vr transition induced by the frequency and strength of the cr_ 
polarized laser field. This dynamic transition is not affected by the strong Coulomb interaction. 
Especially, the spin current appears and is an even function of the phase difference between two 
superconductors when the frequency of the polarized laser field is larger than two times supercon- 
ducting energy gap. Our theory serves as an extension to non-superconducting spintronics. 

PACS numbers: 72.25.Fe, 74.45. +c, 78.67.Hc 
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The interplay between superconducting order and ferromagnetic order has attracted 
considerable attention recently in the proximity effects near the interface of ferromag- 
net /superconductor heterostructures, since it subsumes many fascinating physical phe- 



ications that may complement non- 
and Refs. therein). A interest- 



nomena at the interface and leads to potential af 
superconducting spintronic devices (see Reviews [l|, 
ing physical phenomenon due to the interplay of magnetic and superconducting orders 
is the so-called — tt transition in superconductor-ferromagnet-superconductor (S/F/S) 



junctions 



0, ^, ^, 0, y, 9] and the magnetic quantum dots 10, 11|. The existence of the 



n 



TT junction in layered S/F/S systems was first predicted by Bulaevskii et a/. [4| and Buzdin 
et a/.[5i] to occur for certain thicknesses and exchange field energies of the F layer and was 
later confirmed in the experiment 6[]. These investigations are not only of academic inter- 
est, but also of importance for the solid-state implementation of a quantum bit based on a 
superconducting loop with and vr Josephson junctions [12]. 



The present work is also related to the parametric quantum pumpingjisl. which refers 
to the transfer of electrons coherently between two reservoirs at zero bias. In particular, 
when the charge transfer is quantized per cycle, the pumping could be of importance in es- 
tablishing a standard of electric current [l^ . Therefore, charge pumping has attracted much 
attention since the first proposal by Thoulessjl5|. The first experimental demonstration of 
the pumping of charges was reported by Switkes et a/, [l^. They applied sinusoidal gate 
voltages to an open quantum dot to pump the charge current. More recently, this pumping 
mechanism has been extended to different transport contexts, including the generation of 



spin current 



17, a 
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2l\ and entangled pairs |22l. |23|. 



In this Letter, we present a nonequilibrium charge and spin transport theory of a quantum 
dot under the optical micro-cavity. We find that the laser field acting as the pumping 
forces can induce the dynamic — vr transition which is not affected by the strong Coulomb 
interaction. Our theory goes beyond the adiabatic limit, and is valid for the arbitrary 
pumping frequency and temperature. More interestingly, we find the spin current as an 
even function of the phase difference between two superconductors is pumped when the 
laser frequency is larger than 2 A. 

The system we consider is a quantum dot embedded in a high-Q micro-cavity. Two 
superconducting reservoirs are coupled to the dot via tunneling. The Hamiltonian of the 
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present system reads 

H = Ho + H{t), (1) 

in which the time-independent Hq consists of three parts: Hs describes the left and right 
ordinary BCS superconductors with the energy gap function A^j^ and band width W 

Hs= Yl ^kcl^aCkaa+ Yl l^^i^a^-Ha + ^.c]. (2) 
kcra=L,R ka=L,R 

We assume = A exp(i<l>Q,), i.e., two superconductors have the same energy gap but the 
different phase. stands for the Hamiltonian of single energy level quantum dot with the 
on-site Coulomb interaction 

Hd = 'Y ("(^dld^ + Unini. (3) 

a 

Note that the Zeeman splitting due to the external magnetic field B between two spin 
states is h = ~ = gfisB, and U is on-site Coulomb interaction. Hx is the tunneling 
Hamiltonian between the superconductors and the quantum dot 

^T = Y1 [Tkaacl^Ja + T^^^^Cfeaa] ■ (4) 
kaa 

The most important term H{t) denotes transitions between the different spin states of the 



24| . If we treat the strong laser 



dot, which can be induced by a two-photo Raman process 
field classically, H{t) has the following form 

H{t)=r [d\di exp{-iujt) + h.c.]. (5) 

Here r and uj are the classical Rabi frequency and the cr_ polarized laser frequency, respec- 
tively. This cr_ polarized laser acts as the two pumping forces with the phase difference 
tt/2. It must be pointed out that neglecting the on-site Coulomb interaction and replacing 
the superconductor by normal metal in the above model have been used to generate spin 



current 
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24| . Due to the existence of spin-flip term, we have to work in the general- 
ized Nambu(8)spin space. The standard Keldysh Green's function gives the charge and spin 
current {e = h = 1) 

-^-^Tr{a,,s[G^{E,,E,)J:<{E2) + G<{E,,E^)^l{E^) + h.c]} 
Zn Zn 

*exp[-t{E,- E2)t] (6) 




where (Tc = | | " | , cTg = | | , | . G'^{Ei,E2) and G<(£^i,£'2) are the Fourier 

representations of the retarded and lesser Green's functions in the generahzed Nambu(8)spin 
space with the creation operator = (rfj, d^, rfj, d-^). They have the usual form G^(ti, = 
-i9{ti - t2)({*(ti),*"^(t2)}) and G'<(ti,t2) = i(^"f(t2)^(ti)),rcspcctivcly. S^'<(S) are the 
retarded and lesser self-energies due to the superconducting lead a. The lesser Green's func- 
tion is related to the retarded and advanced Green's functions via Keldysh equation G^ — 
G"'E^G", therefore, the whole problem is reduced to the calculation of the retarded Green's 
function G"". This can be done by two steps: we first calculate the time-independent retarded 
Green's function Gq{E) of the quantum dot, then the full Green's function can be obtained 
exactly from the Dyson equation G^{ti,t2) = GQ{ti — ^2) + / dtGQ{ti —t)H{t)G^{t, 12), which 
has the following Fourier form 

( G;,iE,, E2) G{,{E,, E2)) = [ C2^ Cu ) A/DetA, (7) 

( GUEi, E2) GUEi, E2)) = [ C4. C3. ) B/DetB, (8) 
where the matrices A, B and the coefficients Cij are defined as 

Gl2i{Ei) -Gl22{Ei) j \ Gl^^{Ei - u) Gl^{Ei - u) 




Gij 



(11) 



GS43(^i) -GluiEi) I [ Gl,,{E, + u) Gl,,{E, + u) 

2nGl,^{E,)6{Ei - E2) 
if hi = 1,2 or- ij = 3,4 

2nr[GUEi)Gl,j{E^ - u;) - Gl,2{E,)Gl,^{E, - uME, - E2 - u) 
if i — 1,2 and j — 3,4 

2nr[GUEi)Gl,j{E, + u) - Gl^,{E,)Gl^.{E, + uME, -E2+U;) 
if i — 3,4 and j — 1,2 
Substituting the retarded Green's function into Eq.(6), we can obtain the charge and spin 
current formulae directly, which have the following form 

Ic = J ^IV{<7e[F(£;)/ + - /) + F-{E){f- -f) + h.c.]}, (12) 



dE 
2^ 



Tr{a.[F+(E)(/+ - /) + F-{E){f- - /) + h.c.]}. 



(13) 



Here F{E) = G^{E, E) - S^JE)] + / ^{G^{E, X) [S'^(X) - E^(X)]G«(X, 

= j'^{G^{E,X)T,[T.-{X) - ^^{X)]nG%X,E)T,T.l{E)T^}, F-{E) = 



/|f{G"(E,X)r2[S"(X) - T.'\X)]t2G\X,E)tJ:1{E)ti}^. The matrices n and 











are defined as ri = 






lOi ) 




lOOJ 







respectively. / is the well-known Fermi distribu- 
tion function / = l/(exp(/5£') + l) and = f{E±u!). Eq.(12) shows that the charge current 
comes from two parts: the dissipationless supercurrent Igc = J Tr{o"c[-F(-E)/ + h.c.]} 
involving in the transfer of coherent Cooper pairs and dissipative quasi-particle current 
= / f Tr{ae[F+(E)(/+ - /) + F~{E){r - /) + h.c.]}. Since the Cooper pair is 
singlet, only can the spin current be carried by the dissipative quasi-particle process. We 
must stress that the above charge and spin current formulae go beyond the adiabatic 
approximation, and are valid for the arbitrary temperature and frequency and strength of 
the laser field since we can calculate the retarded Green's function exactly. 

We first consider the non- interacting case on the quantum dot. The time- independent 
retarded Green's function for the quantum dot is given by 



Gl,{E) = {Ei-Hn-J2^';,{E)}-\ 



(14) 



where Hd 



/ et \ 
-e^ 
e| 

\ -et / 



( E -A^ \ 

-a; E 
E Ao 

V a; E / 



with = 1 



for E > —A and = otherwise. Fq, is the linewidth function which is assumed to be a 
constant under wide-band limit. Once having the above time- independent retarded Green's 
function, we can obtain the full Green's function and then charge and spin current by using 
Eqs.(7-13). Fig. 1(a) shows the Josephson current J^c versus the two superconducting phase 
difference $ = — $l of two superconductors for the various polarized laser frequencies lj. 
Our system is in the vr-state when the laser frequency is lower. As the frequency u) increases, 
the TT state will turn into a 0-state first, i.e., a dynamic tt — transition happens. For a 
higher frequency the 0-state will experience a dynamic — vr transition and become a vr-state 
again. This result can be seen clearly in Fig.2, which exhibits the Josephson current Igc as 
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a function of frequency lo. In Fig. 1(b) we plot the Josephson current Igc versus the phase 
difference $ of two superconductors for the different polarized laser strength r. The similar 
0— TT transition can be modulated by the strength r of laser field. Fig. 3 describes the pumped 
quasi-particle charge current Jgc versus the laser frequency. We find that this pumped quasi- 
particle charge current is very small for the lower frequency and increases with frequency 
non-monotonously. The pumped quasi-particle charge current can be understood from the 
co-tunneling process. For example, a spin down electron 1 in the left superconductor tunnels 
into the quantum dot, then a Cooper pair splits into two electrons 2 and 2' with the opposite 
spin. The spin down electron 2' fills in the original position of the electron 1, while spin 
up electron 2 will tunnel into the quantum dot. Finally, the spin down electron 1 and spin 
up electron 2 in the quantum dot change their spin by absorption and emitting a photon, 
respectively, and form a pair to enter the right superconductor. This co-tunneling process 
is depicted in the inset of Fig. 3. Our numerical calculation shows that the pumped quasi- 
particle charge current always satisfies Iqcw — Iqc22 — Iqc33 — Iqcn and gives no spin current 
if only the laser frequency is lower than 2A. However, the pumped spin current appears 
once the laser frequency is higher than 2A. The result can be seen from Fig. 4. Obviously, 
this spin current arises from the additional spin-flip process which is plotted in the inset of 
Fig.4: A spin-down electron from both superconducting leads can tunnel into the quantum 
dot, and due to the (T_ polarized laser field it absorbs a photon and becomes a spin-up 
electron. If the laser frequency is higher than 2A, this spin-up electron can tunnel out of 
the scattering region and goes into the superconducting leads as a quasi-particle. In the 
Fig. 1(c) wc plot the spin current Ig versus the two superconducting phase difference In 
contrast to the Josephson supercurrent, the pumped spin current Is is an even function of 

Finally, we briefly consider there exists a strong Coulomb interacting case {U oo) in 
the quantum dot. We focus on the strong coupling hmit A << Tk, where the Kondo effect 
becomes important and leads to a positive pairing correlation on the quantum dot. This 
situation is well described in the slave-boson language. The physical electron operator 
can be replaced by fa , where b and being the standard boson and fcrmion annihilation 
operators standing for the empty (n| = 0,nj^ = 0) and singly occupied (n| = l,nj^ = 0) 
or (n-f — 0,ni — 1). We then introduce the constriction which prevents double occupancy 
in the quantum dot by means of Lagrange multiplier A. The resulting model can be solved 
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within the mean-field approach. Two constants b and A can be calculated from the following 
self-consistent equations 



Having constants b, X and using Eqs.(7-13), we can obtain the full Green's function and fi- 
nally give the charge and spin current in the strong coupling limit. Fig. 5 shows the Josephson 
current versus the two superconducting phase difference $ of two superconductors for the 
various polarized laser frequencies cu. It is found that as the frequency co increases, the sys- 
tem will go from a 7r-state to a 0-state. This means that the strong Coulomb interaction 
will not affect the dynamic — tt transition. 

In summary, we have presented a charge and spin pumping theory of quantum dot coupled 
to two superconducting leads. The (7_ polarized laser field serving as the pumping forces 
can induce the dynamic — tt transition and spin current. In contrast to the Josephson 
supercurrent, the pumped spin current is an even function of the phase difference between 
two superconductors. The strong Coulomb interaction in the quantum dot does not change 
the dynamic — tt transition behavior. Since the system under investigation is within the 
reach of the present nano-technology, we hope that the present theory can stimulate the 
further study of cavity superconducting electronics and spintronics. 
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Figure Captions 

Fig.l. (a) Josephson current-phase relation for various laser frequencies u. u = 0, 0.5, 
1.5, 2.5 correspond to solid, dash, dash-dot, and short-dash lines, e-f = = 0, r = 0.6, 
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— Fr — 0.1. (b) Josephson current-phase relation for various laser strength r. r — 0.1, 
0.3, 0.5 correspond to solid, dash, dot lines, e-f = —0.2, ej^ = —0.4, u — 0.3, Vl — Vr — 0.1. 
(c) The pumped spin current versus the phase difference $ between two superconctors. We 
have set uj = 2.5 and other parameters are the same as those in Fig. 1(a). Note that all 
parameters are in unit of superconducting energy gap A. 

Fig. 2. Josephson current Isc versus frequency cu. Here ^ — tt/2 and other parameters are 
the same as in Fig. 1(a). 

Fig. 3. Pumped quasi-particle charge current Igc as a function of cu. Here e-f = 0.9, 
ej — 0.5, r = 0.6 and other parameters are the same as in Fig.2(a). Inserted: co-tunneling 
process for the pumped quasi-particle charge current. 

Fig. 4. Pumped quasi-particle spin current Ig as a function of a;. Parameters are the same 
as in Fig. 3. Inserted: spin-flip process for the pumped spin current. 

Fig. 5. Josephson current-phase relation for various laser frequencies cu with a strong 
Coulomb interaction in the quantum dot. u; — 0.6, 1.2, 2.2, 2.8 correspond to solid, dash, 
dot, dash-dot hne. Other parameters are e = —2, = = 0.5, r = 2, $ = 7r/2, 
A = 0.002, W = 100, f3 = 500. 
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